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Quantum entanglement [U [2] plays a vital role 
in many quantum information and communica- 
tion tasks [3j. Entangled states of higher di- 
mensional systems are of great interest due to 
the extended possibilities they provide. For ex- 
ample, they allow the realisation of new types 
of quantum information schemes that can offer 
higher information-density coding and greater 
resilience to errors than can be achieved with en- 
tangled two-dimensional systems (see [4j and ref- 
erences therein). Closing the detection loophole 
in Bell test experiments is also more experimen- 
tally feasible when higher dimensional entangled 
systems are used [5j. We have measured previ- 
ously untested correlations between two photons 
to experimentally demonstrate high-dimensional 
entangled states. We obtain violations of Bell- 
type inequalities generalised to dimensional 
systems [Ij with up to d = 12. Furthermore, the 
violations are strong enough to indicate genuine 
11-dimensional entanglement. Our experiments 
use photons entangled in orbital angular mo- 
mentum (OAM) [7J, generated through sponta- 
neous parametric down-conversion (SPDC) [8, 9j, 
and manipulated using computer controlled holo- 
grams. 

Quantum information tasks requiring high dimensional 
bipartite entanglement include teleportation using qu- 
dits [To", IT , generalised dense coding (i.e., with pairs 
of entangled d-level systems) [12 , and some quantum 
key distribution protocols [13]. More generally, schemes 
like quantum secret sharing [TT , and measurement based 
quantum computation [15 , apply multi-particle entan- 
glement. These are promising applications, especially 
in view of recent progress in the development of quan- 
tum repeaters (see [16j and references therein). How- 
ever, practical applications of such protocols are only 
conceivable when it is possible to experimentally prepare, 
and moreover, detect high-dimensional entangled states. 
Therefore, the ability to verify high-dimensional entan- 
glement between physical qudits is of crucial importance. 
Indeed, much progress has generally been made on the 
generation and detection of high-dimensional entangled 
states (please see [17 and references within). 

In this Letter, we report the experimental investiga- 
tion of high-dimensional, two-photon entangled states. 
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We focus on photon OAM entangled states generated by 
SPDC, and demonstrate genuine high-dimensional entan- 
glement using violations of generalised Bell-type inequal- 
ities [1 . Previously, qutrit Bell- type tests have been 
performed using photon OAM to verify 3-dimensional 
entanglement (see [18 and references within). In ad- 
dition to testing whether correlations in nature can be 
explained by local realist theories [19 , the violation of 
Bell-type inequalities may be used to demonstrate the 
presence of entanglement. Bell- type experiments have 
been performed using two-dimensional subspaces of the 
OAM state space of photons [20l[2T] and experimentalists 
have demonstrated 2-dimensional entanglement using up 
to twenty different 2-dimensional subspaces [22 . Care- 
ful studies have also been carried out to describe how 
specific detector characteristics bound the dimensional- 
ity of the measured OAM states in photons generated by 
SPDC using Shannon dimensionality [23\ 

Our experimental study of high-dimensional entangle- 
ment is based on the theoretical work of Collins et al. [l] , 
which was applied in experiments for qutrits encoded in 
the OAM states of photons [I8l [24] . We encode qudits 
using the OAM states of photons, with eigenstates de- 
fined by the azimuthal index I. These states arise from 
the solution of the paraxial wave equation in its cylin- 
drical co-ordinate representation, and are the Laguerre- 
Gaussian modes I/Gp,^, so called because they are light 
beams having a Laguerre- Gaussian amplitude distribu- 
tion. 

In our setup (Fig.jl]), OAM entangled photons are gen- 
erated through a frequency degenerate type-I SPDC pro- 
cess, and the OAM state is manipulated with computer 
controlled spatial light modulators (SLM) acting as re- 
configurable holograms. Conservation of angular momen- 
tum ensures that if the signal photon is in the mode spec- 
ified by 1^), the corresponding idler photon can only be 
in the mode \—t). Assuming that angular momentum is 
conserved [9 , a pure state of the two photon field pro- 
duced will have the form 

£=00 

where subscripts A and B label the signal and idler pho- 
tons respectively, is the probability to create a pho- 
ton pair with OAM and \t) is the OAM eigenmode 
with mode number I. 

Collins et al. [1 showed that, for correlations which 
can be described by theories based on local realism [1] , a 
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FIG. 1. Schematic of experimental setup for violations of Bell- type inequalities. C{Aa — v^Bh — w) or 

c{e\,e^B) is the 

coincidence count rate when SLM A is in state \v)a or \6\) and SLM B is in state \w)^ or |^^) respectively. 



family of Bell-type parameters Sd satisfy the inequalities 



^(local realism) ^ 



for all d>2. 
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Alternatively, if quantum mechanics is assumed to hold, 
then the violation of an inequality of type Q indicates 
the presence of entanglement. Sd can be expressed as 
the expectation value of a quantum mechanical observ- 
able, which we denote as Sd- The expressions for Sd^ Sd 
and the operators ^2, S^ are provided in Sec. SI of the 
Supplementary Information. 

The parameters 5'^^ are calculated using coincidence 
probabilities for measurements made locally by two ob- 
servers, Alice and Bob, on their respective subsystems, 
which in our case are the signal and idler photons from 
the SPDC source. Alice's detector has two settings la- 
belled by a G {0, 1} with d outcomes for each setting, 
and similarly for Bob's detector with settings h G {0, 1}. 
The measurement bases corresponding to the detector 
settings of Alice and Bob are defined as 
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where v and w both run from to d — 1 and denote the 
outcomes of Alice's and Bob's measurements respectively, 
and the parameters = 0, ai = 1/2, /3o = 1/4, and 
/3i = -1/4. 

The measurement bases ^ind have been 

shown [5l [25] to maximise the violations of inequality ([2| 



for the maximally entangled state of two d-dimensional 
systems given by = ^ J]) j=o b')^ ^ \^)b- It turns 
out that we are able to parametrise these (i-dimensional 
measurement basis states with 'mode analyser' angles Oa 
and Ob^ and write them in the form 

£=+[d/2] 

l^)a = Wa) = ^ E ^'^AOm K), and 
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= ^ E exp[z^^^(^)]K), 



(5) 



= + a/2) 27r/d, 
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The function g{t) is defined as 



9{l) = l- 



{d mod 2)u{£), 



(6) 



(7) 



where [x] is the integer part of x, and u(^) is the discrete 
unit step function. 



Fig. SI shows an example of the experimental data 
points for the self normalised coincidence rates as func- 
tion of the relative angle {6a — Ob) using d = 11 (see 
also Fig. SI in the Supplementary Information). For a 
maximally entangled state 
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FIG. 2. Coincidence count rates (self normalised) as a func- 
tions of relative orientation angles between state analysers 
(Oa — Ob)- Eqn. [ojfor a state with maximal 11-dimensional 
entanglement is fitted to the experimental data with the ver- 
tical offset and amplitude left as free parameters. Errors were 
estimated assuming Poisson statistics. 
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FIG. 3. Experimental Bell- type parameter Sd versus number 
of dimensions d. Sd > ^2 violates local realism for any d > 2. 
The plot compares the theoretically predicted violations by 
a maximally entangled state and the local hidden variable 
(LHV) limit with the experiments. Violations are observed 
for up to d = 12. Errors were estimated assuming Poisson 
statistics. 



where h{i) = 1 for all i when d is odd, and h{i ^ 0) = 1^ 
h{0) = 1 when d is even, the coincidence rate of detecting 
one photon in state \6a) and the other in state \6b) is 
proportional to 



C{OA^OB) = \{OA\{OB\m^ ^ 



cos{d{0/ 



^))-i 
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(9) 

The key result of our paper is displayed in Fig.[3j which 
shows a plot of experimental values of parameter Sd as a 
function of the number of dimensions d. The plot com- 
pares theoretically predicted violations for a maximally 
entangled state, the experimental readings and the lo- 
cal hidden variable (LHV) limit. The maximum possible 



violations (shown in Tab. SI of the Supplementary In- 
formation) are slightly larger than the corresponding vio- 
lations produced by a maximally entangled state. Viola- 
tions persist up to as much as d = 12 when entanglement 
concentration [26^ is applied. We find Su = 2.39 ± 0.07 
and 6'i2 = 2.24 ± 0.08, which clearly violate Sd < 2 (see 
also Tab. S4 of the Supplementary Information). In the 
corresponding experiment using LGp/ modes with only 
p = 0, violations are obtained up to d = 11. Without 
entanglement concentration, we observe violations only 
up to (i = 9 (please see Fig. S2 in Sec. SIV of the Supple- 
mentary Information). Above d ^ 11, the strength of the 
signal becomes so low that noise begins to overshadow the 



quantum correlations. In Fig. [STJ the theoretical predic- 
tion in Eqn. [9] for a state with maximal 11-dimensional 
entanglement is fitted to the experimental coincidence 
data obtained using the mode analyser settings defined 
in Eq. [5]for d = 11, with only the vertical offset and am- 
plitude left as free parameters. The observed fringes are 
seen to closely match those theoretically obtained for a 
state with maximal 11-dimensional entanglement. 

The violation of a Bell inequality in d x d dimensions 
directly indicates that the measured state was entangled. 
It remains to determine how many dimensions were in- 
volved in the entanglement. Measuring the coincidence 
probabilities, i.e., of having the joint state \is) (g) 
(Fig. [4|, together with the parameters Sd for different 
can be seen as a partial tomography of the SPDC source 
state. Numerical investigations indicate that a state hav- 
ing the experimentally observed coincidence probabilities 
and parameters 5*2, S'3, . . . , /Su must contain genuine 11- 
dimensional entanglement. In other words, it is not pos- 
sible to obtain the observed levels of violation with a 
state that contains entanglement involving only 10 di- 
mensions or less. Our analysis assumes a special form of 
the states, based on the coincidence measurement results 
shown in Fig. [4] Further details are given in Sec. SII of 
the Supplementary Information. 

Our results hold much promise for applications requir- 
ing entangled qudits in general. As mentioned earlier, 
progress in the development of quantum repeaters (see 
[16] and references therein) would make quantum key dis- 
tribution using high-dimensional entangled states J3] a 
possible application. Conventional quantum communica- 
tion will fail for sufficiently large transmission distances 
because of loss, and quantum repeaters are one possible 
solution to this problem. Although experimental quan- 
tum key distribution has been demonstrated with 0AM 
qutrits |24^, our findings provide experimental evidence 
that such protocols could be implemented using photons 
entangled in 0AM in up to 11 dimensions, resulting in a 
considerable increase in information coding density. 

A possible extension to our work would be to inves- 
tigate the generation of multi-photon, high-dimensional 
0AM entanglement. One can conceive of achieving this 
using a cascade of down-conversion crystals for generat- 
ing multipartite entangled photons, which has been done 
for polarisation entangled photons fT7|. It also appears 
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FIG. 4. Experimental coincidence rates proportional to the 
probability of measuring the state \ls) \U) with IsAi — 
— 5, . . . , +5. The coloured and greyed-out bars depict the 
measurement results with and without the application of Pro- 
crustean filtering respectively. The measurement time was 
20 s for each combination of Is and 



to be within reach to combine the high-dimensional pho- 
ton 0AM entanglement with entanglement in the po- 
larisation and path degrees of freedom, creating even 
larger hyper-entangled states (please see [17 and refer- 
ences within). 

On a more fundamental note, Bell test experiments 
performed to date have one or both of two main loop- 
holes, namely the locality and detection loopholes. How- 
ever, a recent theoretical work reveals that even low di- 
mensional qudits can provide a significant advantage over 
qubits for closing the detection loophole [5 . In fact, 
it was found that as much as 38.2% loss can be tol- 
erated using 4-dimensional entanglement. Our results 
raise interesting possibilities regarding the role higher- 
dimensional entangled qudits could play in closing this 
loophole. We emphasize that neither the detection nor 
the locality loophole has been closed in our experiments, 
because the overall eflftciency of our experimental setup 
is 1 — 2%, and the switching time for our measurement 
devices (SLMs) is of the order of tens of ms. However, 
closing these loopholes was not the immediate goal of our 
experiments. We are instead using the violation of Bell 
inequalities, up to fair sampling assumptions, as a means 
of verifying the presence of high-dimensional entangle- 



ment, within the framework of quantum mechanics. 

In summary, we have been able to experimentally 
demonstrate violations of Bell-type inequalities gener- 
alised to (i-dimensional systems [Ij with up to d = 12, 
enough to indicate genuine 11-dimensional entanglement 
in the orbital angular momentum of signal and idler pho- 
tons in parametric down-conversion. It appears that this 
could be extended to even higher dimensions by using a 
brighter source of entangled photons. 

METHODS 

In our experiments, we use computer controlled SLMs (Hama- 
matsu) operating in reflection mode with a resolution of 600 x 600 
pixels. In the detection, the SLMs are prepared in the states de- 
flned in Eqns. ([sj respectively. An SLM prepared in a given state 
transforms a photon in that state to the Gaussian \£ = 0} mode. 
The reflected photon is then coupled into a single-mode flbre which 
feeds a single photon detector. Since only the \£ = 0} mode couples 
into the flbre, a count in the detector indicates a detection of the 
state in which the SLM was prepared. The hologram generation 
algorithm introduced in ^7] is applied to conflgure the SLMs. 

Fig.[l]shows the schematic diagram of the experimental setup as 
well as examples of SLM settings used where d = 11. For the SPDC, 
we use a pump beam, with i = 0, produced by a frequency tripled, 
mode-locked Nd-YAG laser with an average output power of 150 
mW at 355 nm. The collimated laser beam is normally incident 
on a 3 mm long BBO crystal cut for type-I collinear phase match- 
ing. A 50:50 beam splitter (BS) then separates the co-propagating 
OAM entangled photons probabilistically into the signal and idler 
paths. Spectral fllters with 10 nm bandwidth are used to reduce 
the detection of noise photons. The coincidence resolving time is 
10 ns and an integration time of 20 s is used for the measurements. 

For tests within a d-dimensional subspace and for odd d, we 
choose the modes £ = —(d— l)/2, . . . , 0, . . . , (fi — l)/2 as the com- 
putational basis states \j) in Eqns. [3] and [4] where j = 0, . . . d — 1. 
For even d, we use i = —d/2, . . . , —1, 1, . . . , c^/2, omitting the 
£ = mode. A projection of the SPDC output state onto a d- 
dimensional subspace results in a non-maximally entangled state 
due to the limited spiral bandwidth [4 . To enhance the entangle- 
ment, we use the so-called Procrustean method of entanglement 
concentration 26 . This is generally done by means of a flltering 
technique which equalises the mode amplitudes, thereby probabi- 
listically enhancing the entanglement of the two-photon state [29] . 
This can be achieved by applying local operations to one or both of 
the signal and idler photons. We choose local operations matched 
to the spiral bandwidth measurement for our SPDC source (please 
see Sec. SIV of the Supplementary Information), so as to ob- 
tain a close approximation to a maximally entangled state. The 
method applied in '30] uses lenses for equalising amplitudes in a 
superposition of three OAM modes. We however use alterations 
of the diffraction efficiencies of blazed phase gratings in the SLMs 
to achieve this goal for up to fourteen modes. Fig. [4] contrasts 
the results of coincidence measurements with and without Pro- 
crustean flltering, with the SLMs in the state {|^a} \^b)} where 
'^aAb ^ {—5,..., +5}. The disadvantage of the Procrustean 
method is the associated reduction in the number of detected pho- 
tons (please see Sec. SIII of the Supplementary Information for 
further details). 
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SI. BELL OPERATOR 

The expression for the generahsed Beh-type parameter derived in [1 can be written as 

[d/2]-l 

Sd= (l-^){ + [^(^o = Bo + A:)+P(Bo=:Ai + A:+l) + P(Ai=Bi + /c) + P(Bi=Ao + A;)] 

fc = 

- [p(Ao = Bo-k-l) + P{Bo = Ai-k) + P{Ai = Bi-k-l) + P{Bi = Ao - k - 1)]}. (SI) 

Here, d is the number of dimensions, and Sd is the Beh parameter, denoted as in yy. The measurement outcomes 
Ai, Bi e {0^ . . . — 1}. P{Aa = Bij) denotes the probabihty that the outcome Aa of Ahce's measurement is the 
same the outcome B^j of Bob's measurement, for the respective detector settings a, 6 G {0, 1}, (denoted in [1 as 
a,6 G {1,2}). That is, 

d-1 

P{Aa = B,) = Y^ P{Aa = j, B, = j). (S2) 

In a similar way, 

d-1 

P{Aa = B,^k)=Y,P[Aa= j, B, = {j + k) mod d] (S3) 

j=o 

is the probabihty that the outcome B^j differs from Aa by /c, modulo d. We also have 

d-1 

P{B, = A, + fc) = ^ P [A, = {j + k) mod B, = j] . (S4) 

Each joint probability is obtained from the photon coincidence count rate C{Aa = i^Bjj = j) divided by the sum 
Ct^cl-, b) of all coincidence rates for a given combination of detector settings for Alice and Bob, 
Cria, b) = Et:r=o C{Aa = i', B, = j'). 

The d basis states correspond to orbital angular momentum (GAM) eigenstates, with an GAM of ITi per photon. 
For odd i runs from —{d— l)/2 to {d — l)/2, and for even £ runs from —d/2 to -\-d/2, but without the £ = 
state. These basis states are denoted as where j = 0, . . . , d — 1. For the signal photon state and odd we choose 
j = £ -\- {d — l)/2^ and for even similarly, j = corresponds to £ = d/2^ j = 1 corresponds to £ = d/2 -\- 1^ and so 
on, until j = d — 1 which corresponds to £ = d/2 (^ = is missing for even d). For the idler photon the ordering is 
the opposite, with low values of j corresponding to high values of £. That is, j = —£ -\- {d — l)/2 for odd d, and for 
even d, j = corresponds to £ = d/2 and so on, until j = d — 1 which corresponds to £ = —d/2 (again without the 
£ = state). In short, a state \£) <^ \ — £) = \ j) ^ \ j) = |j, j), with the connection between £ and j as above. 
As mentioned in the main paper, the Bell-type parameter Sd can be written as the expectation value of a Bell 
operator Sd^ that is, 

Sd{p) = Tr{Sdp). 

Let us define operators P{Aa = Bb + k) as 



P{A, = B, + k) = Y, \r)a\{r + k) mod d)f f {(r + k) mod (S5) 



and operators P(P5 = Aa ^ k) as 

d-1 



P{B, =A, + k) = ^\{r + k) mod d)f |r)f f (r|f ((r + k) mod d\, (S6) 



r=0 
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where the measurement basis states 1'^)^ ^ 
defined in Eqns. 2 to 4 in the main paper. 

From Eqn. [ST| the generalised Bell operator can then be written as 



w)^ for measurement settings a and b and v^w = 0^ . . . - 



1 are as 



[d/2]- 



1 - ^) { + [^(^0 = Bo + /c) + P{Bo = Ai + /c + 1) + P{Ai = Bi + k) + P{Bi 



- [p(Ao = Bo-k- 
For example, the Bell operators for d - 



-l) + P(Bo = Ai- 
2 and d = 3 are 



k) + P{Ai = Bi-k-l + P{Bi = 



Ao + /c)] 
-1)]}. 



(S7) 
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|2,0), 



where the basis states are ordered as |0, 0), |0, 1), |0, 2), . . . , |0, d — 1), |1, 0), |1, 1), . 
\d-l,d-l). 

It is straightforward to check that all diagonal elements of Sd are equal to zero in the 0AM basis |j, /c), as all 
diagonal elements of the operators P{Aa = Bi^ -\- k) and the operators P{B}) = Aa -\- k) are equal to 1/d'^. Also, we 
note that all elements Sjk of the Bell operators, written in the 0AM basis, are nonnegative. This in fact holds for 
all Bell operators until d = 14, and we suspect that it holds generally, although it is not necessary for the purposes 
of this paper to prove it for the general case. 

Clearly, the maximal value of the Bell parameter, that is, the maximal Bell violation, is the largest eigenvalue of the 
Bell operator, and this largest violation is obtained for the corresponding eigenstate. Interestingly, this eigenstate is 
not in general equal to a maximally entangled state. Table SI lists the Bell parameter values theoretically obtained 
for the maximally entangled state \i^)= (l/Vd) Xlj=o b'i)' together with the maximal values of the Bell parameter, 
which, as already stated above, are equal to the largest eigenvalue of the respective Bell operator ^. For any p, the 
corresponding Bell violation is a linear combination of the eigenvalues of Sd, since 



Sd = Tr{pSd) = ^Sk{sk\p\sk), 



(S9) 



where and \sk) are the eigenvalues and eigenvectors of Sd, and therefore ^^i^isklplsk) = Trp = 1. The five largest 
eigenvalues of Sn are shown in Table S2, along with the form of the corresponding eigenstates \sk). 



SIX. EVIDENCE FOR HIGH-DIMENSIONAL ENTANGLEMENT 



Violation of a Bell inequality indicates that the measured state was entangled, but does not directly give 
information of how many dimensions were participating in the entanglement. We will now analyse the results of the 
coincidence measurements of the SPDC output, shown in Fig. 4 of the main paper, together with the obtained Bell 
violations for different (i, and show that the high levels of violations of the tested Bell inequalities indicate that the 
SPDC state is indeed a high-dimensional entangled state. The problem of finding a lower bound on the 
dimensionality required to produce certain quantum correlations has been studied by Acin et al. in [3 , where they 
present examples of correlations that require measurements on quantum systems of dimension greater than two for 
their generation. Although, the depth of analysis presented there would require immense computational efforts for 
11 dimensions, additional experimental results and considerations below allow us to simplify our analysis. 
Assuming conservation of 0AM, meaning that j is always the same for signal and idler photons, both the filtered 
and unfiltered SPDC source states can be written in the form 

d-i 

P=Y1 Cjk\jJ){k,k\. (SIO) 

j,k=0 
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As seen from Fig. 4 in the main paper, for high |-^|, there are nonzero coincidence counts also when j for the signal 
and idler photons differ by ±1. This occurs for approximately 8% of the total count rate. Conservation of angular 
momentum in parametric down-conversion is, however, a theoretically and experimentally well-documented 
phenomenon (please see Ref. [9] of the main paper, and references within). One is therefore justified in assuming a 
state of the form in Eq. ( S1Q[ ), and that any coincidence counts indicating unequal j for signal and idler photons are 



due to imperfect mode selectivity in the measurement, which has greater effect for higher 



A state of the form in Eqn. (SIO), which in addition has only at most {d — l)-dimensional entanglement, can be 



written as a mixture of pure states \tlJm) = Zlj^o ^rnj\j^j)^ where for each iV^rn), it holds that amj = for at least 
one j. By grouping terms with different |V^rn) together, we can further write a state with at most {d — 1) -dimensional 
entanglement as a mixture of no more than d mixed states, p = XI^Zq ^nPn, where each pn is a mixture of states 

IV^rn) for which amn = 0. That is, if pn = Yl^~k=Q Cnjk\j^j){k^ then Cnjk = whenever one or both of j, k are 
equal to n. 

By varying the and the Cnjk in 

d-1 d-1 d-1 

P=^rnpn = ^rn ^ Cn,jk\j J) {k, k\, (Sll) 

n=0 n=0 j,k=0 

one may now investigate what level of Bell violations may be obtained for a state with at most {d — l)-dimensional 
entanglement. Clearly, the maximum violation for such a state is attained when it is pure, if there are no constraints 
in the maximisation. When constraints are included, it is possible in general that there is no pure state which 
satisfies all of them, so that the maximum is obtained with a mixed state. For example, three constraints on the 
Bloch vector of a 2-dimensional state may force it to be mixed. For d > 4 it is in fact enough to consider pure states 
pn = |'0n)('0n|- This follows siucc, for any mixed pn, there is a corresponding pure state with the same diagonal 
components Cnjj, and the same Bell violations for 5*2, ^'3, . . . , Sd- To see this, note that in addition to the d — 1 
diagonal components of pn, we constrain the values of d — 1 Bell parameters, and this constitutes in total 2d — 2 
constraints on p^, which has {d — 1)^ — 1 free parameters. These parameters can be thought of as the diagonal 
components of pn, plus the real and imaginary parts of its off-diagonal elements. We can also think of them as the 
components of the generalised Bloch vector for pn. If the constraints leave at least one component of the generalised 
Bloch vector free, then this last component can be chosen large enough for pn to be a pure state. We thus require 
that (d— 1)^ — l>2(i — 2, which is satisfied for d > 4. For d > 5, we can in fact choose the pn not only to be pure 
but also to have only real matrix elements, since then the number of independent off-diagonal matrix elements in pn 
is strictly greater than the number of Bell inequalities, that is, {d — l){d — 2) /2 > d — 1 if d > 5. This considerably 
simplifies the numerical maximisation, which otherwise would be relatively demanding. 

We focus on the case d = 11, since the level of violation for d = 12 is smaller. The corresponding state space is thus 
spanned by the 0AM basis states \j) = |j, /c), with j, /c G {0, . . . , 10}. For d = 11, we now use the 
experimentally obtained diagonal coincidence probabilities Pjj, with j = 0, . . . , 10, and the experimental values of 
5*2, S's, . . . , 5'io as constraints. All these quantities are allowed to vary within the experimental error bars; errors are 

estimated assuming Poisson statistics for the photon counting processes. By numerical maximisation using 

M athematica® , we find that the largest value of Su, that can be obtained with a state of the form in Eq. (Sll), is 



Sll = 2.14. As mentioned in the main paper, we measure Su to be 2.39 ± 0.07 using all p modes. Using only p = 0, 
we find 2.67 ± 0.22. These and the other experimentally obtained Bell parameters are listed in Tables S3 and S4. In 
other words, the largest violation that can be obtained with a state that has at most 10-dimensional 
entanglement is smaller than the experimentally measured violation, with more than three standard deviations for 
all p modes, and with more than two standard deviations for p = 0. One may therefore conclude, with high 
confidence, that the SPDC output state does contain 11-dimensional entanglement. 
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SIII. PROCRUSTEAN FILTERING 

Let us denote the unfiltered state as pi. The Procrustean entanglement concentration we perform can be considered 
to be a two-outcome generahsed measurement (POM or probabihty operator measure) with measurement operators 
0\0i and 0202- The Hermitian operators Oi and O2 satisfy 0\0i + O2O2 = 1 in the relevant Hilbert space. 
Furthermore, the procedure can be performed using only local operations on the signal and idler beams. The filtered 
state is obtained conditioned on outcome 1, and is given by 

If outcome 2 is obtained, then the filtering has failed, and in our case, we can think of the photon as having been 
absorbed, leaving the vacuum state, 02pi02 /^'^[02pi02)] = |0)(0|. 
Using the filtered state pj, the parameter Sd is 



Sd = TripfSd) = Tr 



(S13) 



The Procrustean filter we applied for all the experiments was by means of local operations of the form 

Oi = Oa ^ Ob^ where Oa — Ob = <3 is approximately a diagonal matrix. For (i=ll, has diagonal elements (1.00, 

0.97, 0.94, 0.92, 0.91, 0.90, 0.91, 0.92, 0.93, 0.95, 0.97). 

Since the purpose of our experiment is to demonstrate entanglement, within the framework of quantum mechanics 
(i.e., assuming that we are testing some quantum mechanical state for entanglement), it is important to note that 
the local filtering can produce neither the appearance of entanglement nor the violation of a Bell inequality from a 
separable state. Any separable state can be written as p = '^^PmPm,A ^ Pm,B- After filtering, we have 

Pf (X Oi pOi = ^^Pm{opm,Ao) {dpm,Bo), which is just another separable state. No separable state can produce a 
violation of any of the inequalities 5'^^ < 2. 

SIV. SUPPLEMENTARY EXPERIMENTAL RESULTS 

Spiral bandwidth. The spiral bandwidth is an important factor affecting the amount of entanglement of the 0AM 
correlated photons for a given selected subspace of the 0AM Hilbert space. The square of a Lorentzian function 

f{£,^) = A^/{^^ + f), (S14) 

gives a good fit to the experimentally measured coincidence rates, which are proportional to |qP and presented in 
Fig. 4 of the main paper. The parameter 7 specifies the half- width at half-maximum (HWHM) of /(^, 7), and A is a 
normalisation constant. We identify 7 with the effective quantum spiral bandwidth. For our source, we obtain 7 to 
be 7.58. As mentioned in the main paper, the finiteness of the spiral bandwidth causes a projection of the SPDC 
output state onto a d-dimensional subspace to result in a non-maximally entangled state [4 . 
Sd versus d with LGp/ modes. Table [S3| shows the values plotted in Fig. 3 of the main paper, using all radial 
modes (all p). This means that the 0AM (denoted by £) states of the photons are used irrespective of their radial 
state (denoted by p), which is another degree of freedom of Laguerre Gaussian beams. The corresponding 
theoretical maximal violations are shown in Table SI. 



Sd versus d using LGp £ modes without Procrustean filtering. Figure S2 (a) shows the experimental results 



(also displayed in Table S5 ) obtained without Procrustean filtering. 
Sd versus d with LGp=Q/ modes. The following results were obtained with the radial index p of the detected 
beams restricted to only p = 0. In this case, the state \t) means \p = 0,^). The data in Table [S4| is plotted in 



Figure S2 (b). This figure compares theoretically predicted violations for a maximally entangled state, the 
experimental readings and the local hidden variable (LHV) limit. Again, the corresponding theoretical maximal 
violations are shown in Table SI. Using only modes with radial index p=0, we observe increased error bars and 
violations only up to (i = 11, due to the reduced count rates resulting from the detection of only p = states. 
Coincidence curves obtained using LGp=o,£ modes. As defined in the main paper, the measurement states in 
terms of 'mode analyser' angles 9a and 9b^ are 

1^^) = — exp[i^^^(^)]K), and = ^ ^ exp [i^B^(£)] K). (S15) 
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Relative orientation of analysers (0^- 0^) Relative orientation of analysers (0^- 0^) 

(a) Experinnent (b) Theory 

FIG. SI. Coincidence count rates (self normalised) as a functions of relative orientation angles between state analysers (Oa — Ob)- 
The figures depict (a) experimental coincidence curves and (b) theoretical prediction for maximally entangled states of two 
d-dimensional systems with mode analyser settings defined in Eg. |S15| using d = 2,5, and 11, as examples. Only modes with 
p = are used here. The observed fringes are typical of genuine 2-, 5-, and 11-dimensional entanglement, respectively. 
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FIG. S2. Supplementary experimental results for Sd versus number of dimensions d using (a) all p modes and no procrustean 
filtering and (b) only p = modes with procrustean filtering. Sd > ^2 violates local realism for any d > 2. The plots compares 
the theoretically predicted violations by a maximally entangled state and the local hidden variable (LHV) limit with the 
experiments. Violations are observed using up to d = 11 dimensions for each qudit. Errors were calculated assuming Poisson 
statistics for the photon counting processes. The measurement time was 20 s per point. 



Fig. [STJa) and Fig. Sl^b) suggest reasonable agreement between the experimental coincidence measurements and 
the theoretical predictions for maximally entangled states. 
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No. of Dimensions d 


Violation for t/?) 


Maximum possible violation 


2 


2.8284 


2.8284 


3 


2.8729 


2.9149 


4 


2.8962 


2.9727 


5 


2.9105 


3.0157 


6 


2.9202 


3.0497 


7 


2.9272 


3.0776 


8 


2.9324 


3.1013 


9 


2.9365 


3.1217 


10 


2.9398 


3.1396 


11 


2.9425 


3.1555 


12 


2.9448 


3.1698 


13 


2.9467 


3.1827 


14 


2.9483 


3.1946 



TABLE SI. Violation of the inequality & < 2 for two dimensional qudits, up to d = 14. This Table shows theoretically 
predicted values obtained for the maximally entangled state (depicted in Fig. 3 of the main paper), and the maximum 
possible violation corresponding to the largest eigenvalue of the Bell operator Sd- 



Eigenvalue of Sn 


Eigenstate 


. . . -dimensional entanglement 


3.1555 


\si) 


= Enio Ci,n\n,n) 


11 


2.4107 




ELo^2,nKn+l) 


10 


2.4107 




En=o^3,n|n + l,n> 


10 


1.9709 




En=o^4,n|n,n + 2) 


9 


1.9709 


\S5)^ 


En=o<^5,n|n + 2,n) 
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TABLE S2. Five largest eigenvalues of ♦S'n and the form of the corresponding eigenstate. 



No. of Dimensions d 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


Parameter Sd 


2.79 


2.78 


2.87 


2.73 


2.76 


2.62 


2.56 


2.46 


2.47 


2.39 


2.24 


2.07 


1.89 


Standard deviation a 


0.03 


0.04 


0.04 


0.05 


0.06 


0.07 


0.07 


0.07 


0.07 


0.07 


0.08 


0.08 


0.08 



TABLE S3. Experimental Bell- type parameter Sd as a function of the number of dimensions d. All radial modes (all p) are 
used here, and violations of *S'cz < 2 are observed for up to (i = 12. 



No. of Dimensions d 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


Parameter Sd 


2.45 


2.4 


2.67 


2.46 


2.79 


2.71 


2.65 


2.7 


2.54 


2.67 


2.1 


2.11 


1.69 


Standard deviation a 


0.09 


0.1 


0.11 


0.12 


0.14 


0.14 


0.16 


0.2 


0.21 


0.22 


0.2 


0.22 


0.24 



TABLE S4. Experimental Bell-type parameter Sd as a function of the number of dimensions d. Only 0AM modes with p = 
are used here, and violations of *S'd < 2 are observed for up to d = 11. 



No. of Dimensions d 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


Parameter Sd 


2.76 


2.77 


2.71 


2.69 


2.53 


2.49 


2.31 


2.19 


1.95 


2.05 


1.75 


1.65 


1.32 


Standard deviation a 


0.03 


0.04 


0.04 


0.05 


0.05 


0.06 


0.06 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 



TABLE S5. Experimental Bell- type parameter Sd as a function of the number of dimensions d without entanglement concen- 
tration. Violations of ^Scz < 2 are observed for up to d = 9. All p modes are used here. 
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